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Abstract 

We derive several entanglement conditions employing non-hermitian 
operators. We start with two conditions that were derived previously 
for field mode operators, and use them to derive conditions that can be 
used to show the existence of field-atom entanglement and entanglement 
between groups of atoms. The original conditions can be strengthened by 
making them invariant under certain sets of local unitary transformations, 
such as Gaussian operations. We then apply these conditions to several 
examples, such as the Dicke model. We conclude with a short discussion 
of how local uncertainty relations with non-hermitian operators can be 
used to derive entanglement conditions. 



1 Introduction 

Entanglement has shown itself to be a valuable resource in quantum information 
in applications ranging from communication protocols, such as dense coding, 
to quantum computing. This has led to a substantial effort to understand and 
characterize entanglement (see [1] and [2] for recent reviews) . There is no simple 
universal test that enables one to tell whether a given state is entangled, but 
there are many sufficient conditions. These include, for example, entanglement 
witnesses. An entanglement witness is an operator whose expectation value is 
nonncgative for separable states, but its expectation value can also be negative, 
and the states for which it is are entangled. For continuous-variable systems, 
there are entanglement criteria involving the expectation values of powers of 
creation and annihilation operators [3]-[S]. Here we would like to expand upon 



the work in [7] to find stronger entanglement conditions, and to find conditions 
for entanglement not just between field modes, but between atoms and field 
modes or between groups of atoms. 

In [7] two conditions that enable one to determine whether the modes in 
a two-mode state are entangled were derived. These conditions are that the 
modes are entangled if either 

\{A ] B)\ 2 > (A^AB^B), (1) 

or 

\(AB)\ 2 > {A*A){BtB), (2) 

where A is any power of the annihilation operator for the first mode and B is 
any power of the annihilation operators for the second mode. These conditions 
are sufficient, but not necessary, to demonstrate entanglement, that is if cither 
one is satisfied, the state is entangled, but if neither one is satisfied, we cannot 
say anything about the entanglement of the state. 

The derivations of these inequalities made no use of the special properties 
of annihilation operators. In fact, if we have a bipartite system described by 
a Hilbert space Ti = TC a ® Tib, and A is any operator on Tt a and B is any 
operator on Hb, then if a state satisfies either of the above conditions, then 
that state is entangled. If A or B is hermitian, these conditions are useless, 
because the Schwarz inequality guarantees that they will be violated for all 
states. Consequently, we will be studying the detection of entanglement with 
non-hermitian operators. 

There have been some previous studies of entanglement making use of non- 
hermitian operators. Toth, Simon and Cirac derived an inequality based on 
uncertainty relations with the number operator and the mode annihilation op- 
erator that can detect entanglement between modes As was mentioned in 
the previous paragraph, Zubairy and one of us derived a set of entanglement 
conditions based on mode creation and annihilation operators, and a very gen- 
eral set of conditions based on the same operators was derived by Shchukin and 
Vogel 0. 

We will also show how our original entanglement conditions can be strength- 
ened by making use of local unitary invariancc. In the next section this will be 
studied for some simple, low-dimesnional systems. We then go on to look at 
an atom, or atoms coupled to a single-mode field and develop entanglement 
conditions for the field-atoms system that are invariant under Gaussian trans- 
formations on the field mode. In the following section, we demonstrate that all 
of the entanglement conditions we find can be derived from the partial transpose 
condition, i.e. that if the partial transpose of a state is not positive, then the 
original state is entangled. However, the conditions we derive are much easier 
to use than the partial transpose condition itself. We then go on to study the 
application of invariance under Gaussian operations to entanglement conditions 
for two modes. We then examine entanglement in two extended examples, the 
Dicke model and light modes coupled by two beam splitters. In both cases we 
find a connection between sub-Poissonian statistics of an input light mode and 
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subsequent entanglement in the system. Finally, we show how a non-hermitian 
version of local uncertainty relations can be used to derive further entanglement 
conditions between field modes, and between a collection of atoms and a field 
mode. 



2 Correlated subspaces 

As mentioned in the Introduction, another clement we will be using in our study 
of entanglement is its invariance under local unitary operations. Because of this 
invariance, it would be nice if our criteria for entanglement also possessed this 
invariance. As we shall see, we cannot always completely satisfy this condition, 
but we can often satisfy it for some subset of local unitary operators. 

Let us begin by considering the first inequality above. Let \ai) a and |a 2 ) a 
be two orthogonal vectors in Tt a and \fli)b and |/3 2 )b be two orthogonal vectors 
in TL b . Now let A = |a 2 )(ai| and B — |/3i)(/3 2 |, and consider the state 

\ip) = ci|ai) |/9i>6 + c 2 |a 2 ) a |/3 2 ) h , (3) 

where |ci| 2 + | C2 1 2 = 1. We find that 

\(A^B)\ 2 = |c lC2 | 2 , (A^AB^B) = 0, (4) 

so that Eq. ([!]) is a good way to detect the entanglement of \tp). 
It will still detect it if we add noise to the state. Defining 

2 2 
Pa=X>iX«il. P,3 = ElWil, (5) 

and considering the density matrix 

p=s\i>){i>\+ ^^Pa^Pp, (6) 

where < s < 1, we find that Eq. shows that the state is entangled if 

s> — — ■ (7) 

If c\ = C2 = l/v2 this gives s > (v5 — l)/2 and if |cic 2 | <C 1, then we have 
s > 1 — 4|cic 2 | 2 . Another possible density matrix that can be detected by this 
condition is one of the form 

p = s\i/j){il;\+p 0l (8) 

where the vectors \otj) a \Pk)b, f° r i> k — 1,2 are in the null space of pa, and 
Tr(po) = 1 — s. For this density matrix, Eq. ([IJ will show that it is entangled if 
s > 0. 
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Now let us go to a more complicated situation. Let Si be the linear span 
of the vectors {|ai), \ct2)}, and 5*2 be the linear span of the vectors {|«3), (0:4)}, 
where the vectors {\aj) G H a \j = 1, ... ,4} form an orthonormal set. We want 
an entanglement condition that will detect entanglement in vectors of the form 

\1>) = ^=(\vi)a\0l)b+\V2)a\fo)b), (9) 

where \vi) E Si and \v 2 ) € S 2 . We would also like the condition to be inde- 
pendent of the specific vectors \vi) a and \v 2 ) a . One way to approach this is to 
make use of the fact that entanglement is invariant under local unitary trans- 
formations. Suppose that U a is a unitary operator on 7i a that leaves Si and 
S 2 invariant. For any two vectors in Si, \vi) a and \v'i) a , and two vectors in 52, 
\v 2 ) a and \v' 2 ) a , we can find a U a such that U a \vi) a = \v[} a and U a \v 2 ) a = W 2 ) a - 
Therefore, if we have an entanglement condition that is invariant under the 
transformations, U a , we will have one that is independent of the vectors \vi) a 
and \v 2 ) a in the equation for \tp) above. We can find such a condition by noting 
that if we have an operator of the form |o!3) (ai| that maps Si to S 2 , then 
U a \a^) a {oii I C^qT 1 wm be a linear combination of the operators \aj) a (a.k\, where 
j = 3, 4 and k — 1,2. Therefore, let us choose 

A = zi I a 3 ) a (ai I + z 2 \ a 4 ) a (a>i\ + z 3 1 a 3 ) a (a 2 1 

+z 4 \a 4 ) a (a 2 \, (10) 

where the complex numbers Zi, . . . , z 4 are arbitrary. The operator B is as before. 
We then have that 

4 

\(A^B)\ 2 - (A^AB^B) = z*jM jk z k , (11) 

where the 4x4 matrix M depends on the state being considered. If M has a 
positive eigenvalue, then by choosing the Zj, j = 1, . . . , 4 to be the components 
of the corresponding eigenvector, we will have an operator A that shows that 
the state we are considering is entangled. 

Let us carry this out explicitly for the density matrix 

p=s\i>){il>\+ 1 ^ S P a ®P , (12) 

where P a is the projection onto Si U S2 and Pp is as before. We then find that 

s 2 1 - s 
M jk = -jVjVl g— Sjk, (13) 

where 

Vi = {vi\ai)(a 3 \v 2 }, i] 2 = (vi\ai) (a 4 \v 2 ) , 

773 = (vi\a 2 ) (a 3 \v 2 ) , 774 = (vi\a 2 )(a 4 \v 2 ). (14) 
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We can then express M as 



M = j\v)(v\- 



1 - s 



(15) 
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where the vector \rj) has components given by rjj, j = 1, . . . , 4. It is then clear 
that M has three negative eigenvalues, corresponding to directions orthogonal 
to \rj), with the remaining eigenvalue equal to (2s 2 + s — l)/8. This is positive 
if s > 1/2, and therefore, the state is entangled if s > 1/2. Note that this 
condition is independent of \vi) a and \v 2 ) a - 

3 Two-level atom coupled to the field 

The above entanglement condition can be applied to a two-level atom coupled 
to a single-mode field. The atom can either absorb a photon and go from its 
lower to its upper state, or emit a photon and go from its upper to its lower 
state. Let the lower state of the atom be \g) and the excited state be |e), and 
let us consider photon states with 0, 1, 2, and 3 photons. If we start the atom 
in its excited state and a superposition of 0, and 2 photons, as the system 
evolves there will be some amplitude to be in these states, but there will also 
be amplitudes to be in the states |<?)|1) and |#)|3). The upper state of the atom 
will be correlated with photon subspace spanned by |1) and |3), and the lower 
state will be correlated with the subspace spanned by |0) and |2). Therefore, 
the entanglement conditions developed in the previous section could be used to 
test for entanglement in this system. 

In considering an atom interacting with a single-mode field, we do not usually 
confine our attention to states of three photons or fewer, so a different set of 
entanglement conditions could be more useful. Instead of limiting the photon 
number, we will consider all possible photon states. In that case we will have 
to give up invariance of the conditions under all unitary transformations of the 
photon Hilbert space. We can, however, derive some simple conditions if we 
restrict the set of unitary transformations to those corresponding to Gaussian 
operations. These operations consist of translations, rotations and squeezing 
transformations. In particular, if a and a* are the annihilation and creation 
operators for the mode, the translations are given by the operators D(a) = 
exp(aa^ — a* a), where a is an arbitrary complex number, rotations are given 
by R(9) — exp(i9a^a), where < 9 < 2ir, and the squeezing transformation 
is given by S(z) — exp[(z*a 2 — z(a^) 2 )/2], where z = rexp(i<f)) is a complex 
number. These transformations act as follows: 



Note that these transformations send the creation and annihilation operators 
into a linear combinations of the annihilation operator, the creation operator, 
and a constant. 



D(a)^aD(a) 
R{9)^aR{9) 
S(z)^aS(z) 



= a + a, 
id 

= e a, 

= acosh r — a) e~' 1 ^ sinh r. 



(16) 
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In order to find an entanglement condition that is invariant under Gaussian 
transformations of the field mode, we set 

A = a-=\g){e\, 

B = zi(a- (a)) + z 2 (a^- {a% (17) 
and substitute these into Eq. ([!]). The result is 

2 

ZjM jk z k > 0, (18) 



where now 
M = 



\(a+Aa)\ 2 - (P e (Aa)tAa) (ff-(Aa)t){ff+(Ao)t) - (P e ((Aa 
a- Aa) (a+Aa) - (P e (Aa) 2 ) \ (a- Aa) \ 2 - (P e Ao(Aa)+) 

(19) 

where Aa = a — (a), cj + = |e)(p| and P e = \e)(e\. If we transform the state of 



the field mode with a Gaussian transformation, the effect in Eq. ( 18 ) is just to 
change the values of z± and z 2 . 



Now if we can find a value of z\ and z 2 so that Eq. ( 18 1 is satisfied, then 
the state is entangled. This will be possible if the matrix M has a positive 
eigenvalue. Therefore, our entanglement condition, which is now invariant under 
Gaussian transformations of the field mode, is that the state is entangled if 
M has a positive eigenvalue. That means, for example, that it can detect 
entanglement in states of the form I at ® S(z)D(a)(\e) |0) + \g)\l))/y/2, where I a t 
is the identity operator for the atomic system, for any value of a and z. 

Let us now look at two examples of atom- field entanglement. Let us first 
consider a two-level atom in the rotating wave approximation interacting with 
a single field mode, which is initially in a thermal state. On resonance, the 
Hamiltonian for this system is 

H = wa)a + ^<t z + n{a + a + cr~a t ). (20) 

The thermal state density matrix is given by 



^ e ™ = — y£ [—) \n)(n\. (21) 



(l + n) 

^ ' 71=0 

With this initial state we find that the off-diagonal elements of M arc initially 
zero, and remain zero for all times. The element in the lower right-hand corner 
(M22) is always negative, so only the clement in the upper left-hand corner 
(M11) can possibly become positive. We plot M u for several different values of 
n in Figure 1. As can be seen, for sufficiently small values of n, our condition 
shows that there are times at which the atom and the field mode are entangled. 

Let us now look at two two-level atoms interacting on resonance with a single 
mode field. The Hamiltonian of the system is now 

H = uoa)a+ —(erf + of) + K {{ a(J t + ^^l) + ( a(J 2 + a 1 ^))- ( 22 ) 
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Figure 1: The eigenvalue Mn as a function of time for average photon numbers 
n = 0.03 (solid line), n = 0.02 (dashed line), and n = 0.01 (dotted line). The 
inset shows Mn as a function of n at Kt = 2.24. 



This Hamiltonian preserves the total number of excitations of the system, so we 
can express the total Hilbert space of the system as a direct sum of subspaces 
containing fixed numbers of excitations, i.e. H = (BTt^ . For n > 1 the 
subspace Ti^ is spanned by the vectors \n,g,g), \n — 1, g,e), \n — l,e,g), and 
\n — 2, e, e), where the first slot contains the photon number, and the remaining 
two slots contain the states of the atoms. For n = 1 it is spanned by \l,g, g), 
\0,g, e), and |0,e,<7). If we start in the state \n,g,g) at time t = 0, we find that 
M12 = M 2 i = for all times, and M 2 2 < 0. Thus our entanglement condition 
reduces to Mn > 0. Finding an explicit expression for Mn (see the Appendix) 
this condition becomes 

- sin 2 (m)[cos(m) + 2(n - l)] 2 

-{n - 1) [2(n - 2)(cos(fit) - l) 2 + (2n - 1) sin 2 (m)j > 0, (23) 

where = Ky/2{2n — 1). It is worthwhile to look at some special cases. For 
n = 1 and n = 2, the condition reads as 

sin 2 (fit) cos 2 (fit) > 0, (24) 

| cos(m) + 2| > 4= i (25) 
v2 

respectively. For the single-photon case, it is not known which of the two atoms 
absorbs the photon, so that the atoms are always entangled. The condition 
becomes more restrictive as n increases. For times in the neighborhood of 2kn/Q 

t = e + k — , (26) 
where e is small but nonvanishing and k is an integer, one can make an e- 
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expansion of the left-hand side of Eq. ( 23 1 around to obtain 



(2n - l)e 2 - -(3n 2 + n + 4)e 4 + 0(e 6 ) > 0. 



(27) 



This can be satisfied for an arbitrarily large n, provided e is small enough, that 
is there always exists a time interval about 2kir/£l during which entanglement 
is detected. Obviously, for an increasing n, this time interval is reduced. 

We can also study the entanglement between both atoms and the field. Now 



we choose A = a and B = J with J 
for entanglement becomes 



<J l + a 2 in Eq. (Il| . The condition 



\(aU-)\ 2 - (aW+J") > 0, 
or, by using the state at time t as given in the Appendix, 
sin 2 (m)[cos(f2i) + 2(n - l)] 2 



(28) 



2n- 1 



-(n- 1) (ra-2)(cos(fit) - l) 2 + (2n - 1) sin 2 (fit) > 0. (29) 



For n = 1 and n = 2 the above equation is the same as Eq. ( 23 1 . In general 



this condition is easier to satisfy than the condition ( 23 1 due to the fact that 



its second term is smaller than its counterpart in Eq. ( 23 1 . As a consequence 



there are moments when we detect entanglement between the field and both 
atoms, while no entanglement is detected between the field and a single atom. 



For times about 2&7r/f2, see Eq. (26 1, one can expand (29 1 with respect to e 
around 

(2n - l)e 2 - ^(3n 2 + Hn + 2)e 4 + 0(e 6 ) > 0. (30) 

Again for large n one can detect entanglement during time windows centered 
around 2k'ir/Q. 



4 Relation to partial transpose condition 

The derivations of the entanglement conditions, Eqs. ([IJ and ([2]), which were 
presented in [TJ, did not make use of the Positive Partial Transpose (PPT) 
condition. This condition states that the partial transpose of a separable density 
matrix is a positive operator, which implies that if the partial transpose of 
a density matrix is not positive, the original density matrix was entangled. 
However, it was subsequently shown that in the case that A and B are powers 
of mode creation and annihilation operators, the conditions are a consequence 
of the PPT condition [ID]. We will now show that this is true in general, i.e. 
that no matter what the choice of A and B, the entanglement conditions in Eqs. 
and ([2]) are consequences of the PPT condition. 

Let us begin by discussing one of these conditions. Our system is divided 
into subsystems a and b, and let A be an operator acting on subsystem a and B 
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an operator acting on subsystem b. A separable state must satisfy the condition 

m 

\{A^B)\ 2 < (A^AB^B). (31) 

If this condition is violated, the state is entangled. Now let us see if we can 
derive this condition from the PPT condition. We start with the condition 

\(AB)\ 2 < (A^AB^B), (32) 

which follows from the Schwarz inequality. These expectation values are taken 
with respect to a density matrix p, and the matrix element of this density matrix 
are 

= ( a {m\ b {(i\)p(\n) a \v) b ), (33) 

where {|™) a } is an orthonormal basis for a and {|/i)t>} is an orthonormal basis 
for b. Now define the partial transpose of p with respect to subsystem a, p Ta , 
to be the operator with matrix elements 

Pm,fi;n,u ~ Pn,fi;m,v- (34) 

For a separable density matrix, p Ta should also be a valid density matrix, that 
is it should be positive and have a trace equal to one. That implies that the 



inequality in Eq. ( 32 1 will hold if the expectation values are taken with respect 
to p Ta . We have that 

Tr(ABp Ta ) = Y,Y, A ™ B »A;nv 

m,n 

= J^(^-Ln) Bvfj,Pn,iJ,;m,v- (35) 

m,n fj,,i/ 

If A\ nn is real, then we see that Tr(ABp Ta ) = Tr(A^Bp). This would be the 
case, for example, if A were a product of mode creation and annihilation opera- 
tors and the basis is the number-state basis. Similarly, we find that if (A' A) mn 
is real, then Tr( J 4 t AB t Bp T ») = Tr(A^AB^Bp). If the density matrix is sepa- 
rable and all of the reality conditions are satisfied, then we have that Eq. d32J, 



with the expectation values taken with respect to p a , implies Eq. (311, with 



the expectation values taken with respect to p. Under these conditions, the 



inequality in Eq. (31 1 follows from the PPT condition. 

As we can see, the above derivation depends on the fact that certain matrix 
elements are real. The condition can be relaxed a bit; as long as the matrix 
elements A\ nn all have the same phase, the derivation will work. However, the 



original derivation of the Eq. ( 31 1 did not require these conditions on the matrix 
elements. 



It is the none the less the case that Eq. (31 1 can be derived from the PPT 



condition with no additional assumptions. Consider two operators F acting on 
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H a and B acting on Tib- We have that, as in the previous paragraph, for a 
separable state 

\(FB}\ 2 < (F^FB^B). (36) 
Now, from above, we see that 

Tr(FBp Ta ) = Y,H^nYB^ Pnm , (37) 

m,n fijV 

Tv(F^ FB^ Bp Ta ) = ^^[(FtF) mn ]*(i?tB)^p n ^ ;m ^. (38) 

m,n \x,v 

Define the operator G, acting on 7Y a , by 

G = J2KJm)(n\. (39) 

m,n 

The quantities involving the partially transposed density matrix can then be 
expressed as 

Tr(FBp T °) = {&B), (40) 
Ti(F^FB^Bp T -) = (G^GB^B). (41) 

For any separable density matrix, we must have 

\(G*B)\ 2 < (G^GB^B). (42) 
In order to recover our original inequality, choose 

F = Y,{A mn y\m){n\, (43) 

m,n 

which implies that G — A, and this completes the derivation of the inequality 
from the PPT condition. The derivation, for separable states, of the condition 

|(AB}| 2 < (A^A)(B^B), (44) 

using the PPT condition, is similar. 



5 Two modes 

Let us now look at the case of two modes with annihilation operators a and b. 
We will first find an entanglement condition that is invariant under Gaussian 
transformations of one of the modes. Let us set A = zi(Aa)^ + z 2 Ao and B = b. 
Then the condition \(A^B)\ 2 > (A^AB^B) can be written as 

(v\M\v) > 0, (45) 

where \v) is the two-component vector 

!«>=(*)' ( 46 ) 
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and M is the 2x2 matrix 

_ / |<Aa6)| 2 - (Ao(Ao)t6t6) (Aab)({Aa)^b)* - ((Aa) 2 tfb) \ 
~ \ ((Aa)^b)(Aab)* - ((Aa^tfb) |<(Aa+6>| 2 - ((Aa)tAa&t&) J ' 

(47) 

In this form, we can see that if M has at least one positive eigenvalue, we can 
find a vector \v) so that the entanglement condition (v\M\v) > is satisfied, so 
that we can say that a state is entangled if the matrix M that results from it 
has a positive eigenvalue. As before, this condition is invariant under Gaussian 
operations, and, consequently, so is our entanglement condition. For this simple 
2x2 case, we can make the condition more explicit. Denoting the matrix 
elements of M by Mjf., where j, k = 1,2, the larger of the two eigenvalues is 

A = \{{M n + M 22 ) + [(M la - M 22 ) 2 + 4|M 12 | 2 ] 1/2 }, (48) 

and the condition that A > is either that (Mn + A/22) > 0, or that |Mi 2 | 2 > 
M n M 22 . 

As a short example of how our new condition is stronger than the one orig- 
inally proved in [TJ, consider the state 

W) = ^S a {z) <g> I 6 (|0>.|1> 6 + |l)«|0) t ). (49) 

This state is obtained by applying a Gaussian operation to the state — 
(|0) a |l)& + |l)a|0)&)/V2 . For the state \ip), the matrix in the previous para- 
graph is diagonal, and the lower right-hand element is positive, so that the state 
is entangled. Since \ip') is obtained from \ip) by a Gaussian operation, the cri- 
terion in the previous paragraph will also show that \ip') is entangled for any z. 
However, if we apply the criterion 

\(a*b)\ 2 > (atatfb), (50) 

presented in [7J, we find that it only shows that the state is entangled if tanh \ z\ < 
1/V2. Therefore, the new condition is an improvement on the old one. 
Now suppose we set 

A = zx{a - (a)) + z 2 (a t - (a 1 )), 

B = w 1 {b-(b)) + w 2 (tf-(tf)). (51) 
Let Hi and H2 be two two-dimensional Hilbcrt spaces with 

eH u (52) 



e n 2 . (53) 



and 



Z2 



\v) 



Wl 

w 2 



11 



The condition \(A^B)\ 2 > (A^ AB^ B) can be written as ((u| <g> (v\)X(\u) <g) \v)) > 
0, where X is a linear hermitian operator on Tix^Tiz, and this operator depends 
on the state being considered. Therefore, a state is entangled if there exists a 
product state, \u) ® |u) in Hi <8> H2, such that 

((u| <8> («|)JT(|u) ® > 0. (54) 

One way of determining whether the above condition can be satisfied is by 
examining the operator X v — (v\X\v), where \v) is a vector in 7i 2 and X v is 
an operator in Hi. If for any vector \v) £ Hi the operator X v has a positive 



eigenvalue, then the condition in Eq. ( 54 1 can be satisfied. This can be somewhat 



tedious since we have to consider all possible vectors \v) , so having some simpler 
criteria is also useful. 

The first criterion involves the eigenvalues of the reduced matrix Xi = 
Tr 2 (X). We will show that if Xi has a positive eigenvalue, then a product 
vector satisfying the above equation exists. Suppose that Xi has a positive 
eigenvalue, A, and that the corresponding eigenstate is \u). In addition, let 
{l w j)b = 1 , 2} be an orthonormal basis for H.2- We then have that 

2 

(u\Xi\u) = 5^((t*| ® {vj\)X(\u) <g> |^>) = A > 0. (55) 
3=1 

Since the sum in the above equation is positive, at least one of the terms must 
be positive, i.e. ((u\ ® (vj\)X(\u) ® \vj)) > for some j. Therefore, if Xi has 
a positive eigenvalue, the state is entangled. This is, however, only a sufficient 



condition for a product vector satisfying Eq. (54 1 to exist, not a necessary one 



That is, it is possible for Xi to have only negative or zero eigenvalues and still 



be able to find a product vector satisfying Eq. (54 1. An example is provided by 
the operator 

X = (|01> + |10»«01| + (10|) - 2(|01) - |10»«01| - (10|), (56) 

for which Xi = -I, but ((+x\(+x\)X(\+x)\+x)) > 0. 

A second criterion involves the eigenvalues of X itself. If X has two or more 



positive eigenvalues, then we can find a product vector satisfying Eq. (54). This 
can be seen by showing that a product vector can be found in the subspace 
spanned by the eigenvectors corresponding to the two positive eigenvalues. Let 
\xi) and I22) be the two eigenvectors with positive eigenvalues, and let the 
Schmidt decomposition of |xi) be given by 

2 
i=i 

We can expand \x2) in the Schmidt basis of \xi) as 

2 2 

3=1 *!=i 



12 



The vector a|xi) + f3\x 2 } is then given by 

2 2 

3=1 fc=l 

where 

Cjk = aKiSjiSki + an 2 6 j2 5 k2 + Pd ]k . (60) 

Now, a\x\) + f3\x 2 ) will be a product state if [cxi/^vz) = (C21/C22) , and making 
use of the above equation, this condition can be expressed as 

yKx+dxx d 21 , . 
w = ' ( 61 ) 

«12 2/K 2 + d 2 2 

where y = a/ (3. This leads to a quadratic equation for y, which can be solved. 
Therefore, there is a product vector in the span of |xi) and \x 2 ). 

A simple example of this procedure is given by examining the entanglement 
of the state 

P = <#oi)<V>oi| + ® Poi\ (62) 

where |V>oi) = (|0) a |l)& + |l) a |0) 6 )/\/2, and P ( f and P ( ^ project onto the 
zero and one-photon states of the a and b modes, respectively, i.e. P^f* = 
|0) a (0| + 1 1) a, (1 1 , and similarly for '. In Ref. [7] it was found that using the 
condition in Eq. ([!]) with the choice A — a and B — b shows that this state is 
entangled for s > (y5 — l)/2. With the choice of A and B given above, we find 
that the eigenvalues of X(v) are given by the solutions of the equation 

A " ~ (^) A + ~ S3(K|2 " |W2|2)2 ~ 2fi2(K|4 + |W2|4)] = °' (63) 

where we have set |i«i| 2 + |w 2 | 2 — lj because only the direction of \v) is impor- 
tant. One of the eigenvalues will be positive if the last term in this equation is 
negative, and this gives us the condition 

1 < s 3 (K| 2 - K| 2 ) 2 + 2s 2 (K| 4 + \w 2 \ 4 ). (64) 

From this we find that the state is entangled if s > 0.474, which is an improve- 
ment over (a/5 - l)/2 = 0.618. 

6 Extended examples 
6.1 The Dicke model 

We would now like to use the simplest versions of our entanglement conditions 
to study entanglement between different subsystems in the Dicke model. The 
entanglement of the ground state in the Dicke model was studied in [11], but 
we will concentrate on entanglement generated by the dynamics. This model 
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consists of N two-level atoms, which are contained in a volume whose dimensions 
are small compared to an optical wavelength, coupled to a single mode of the 
radiation field. Assuming again that the interaction is on resonance, and that 
there is no atom-atom interaction, the Hamiltonian of the system can be written 
as 

N N 
i=l i=l 

= uj^a + iuS z + k(ciS + + a^S~), (65) 
after the introduction of the collective spin operators 

JV N 

4=1 1=1 

S* = ^[S+S-]. (66) 

The equations of motion resulting from this Hamiltonian arc difficult to solve, so 
we will introduce an approximation based on the Holstein-Primakoff represen- 
tation of the spin operators. The Holstein-Primakoff transformation expresses 
the collective spin operators in terms of the bosonic operators £t and £. If we 
choose the ground state of these new operators to correspond to the atomic state 
in which all of the atoms are in their lower energy level, the transformation is 
given by 

S + = ZW2S-eZ, S- = ^2S - tftt (67) 

with S = N/2. From these definitions, it follows that 

S z = -S + &£. (68) 

If the number of excitations of the atomic system remains small with respect 
to the total number of atoms, i.e. -C N, then we can expand the square 

roots keeping only the lowest order terms 

S+ w Vn$\ S- w y/N£, (69) 

and the Hamiltonian of the system simplifies to 

H = uja i a + uj£,^ + KVN(a^ + a^)- (70) 

So far, we have grouped the atoms into a single subsystem. However, we 
would like to study the entanglement between groups of atoms, so we will split 
them into two groups, each group having its own collective spin operators. Let 
us divide the atoms in two groups, one consisting of k atoms and the other 
consisting of N — k atoms, 

st =£ti<tf«v^, s+ = Y^ =k+l ^^VN^kil 

Si =£!u*r«VA£i> s 2 " = ££ fc+ i*r«v^fe&- (7i) 
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In the Holstein-Primakoff representation with only the lowest order terms re- 
tained, the Hamiltonian of this tripartite system is 

H = w(o+a + + 46) + nVk{a£,l + a^i) + Ky/N - k(a^ 2 + a%). (72) 

We would now like to diagonalize this Hamiltonian. In order to do so, we 
first express it as 

H = vHlv (73) 

with v — (a,£,i,^2) T and 

(u> Kyk k\/ N — k \ 

nVk w . (74) 

nVW^k ui J 

The matrix M can be diagonalized by the introduction of the new modes 



1 / k , N-k, 
lN~k^ [k * 

^ = y—^-y^ 

b 2 = —a+\ £i + \ £2, (75) 

to which correspond the eigenvalues 

Ao = to — kvN 
Ai = w, 

A 2 = uj + nVN, (76) 
respectively. Setting f2 = ny/N, the Hamiltonian can now be expressed as 

H=(u- n)blb + ub\bi + (lo + Q)b\b 2 . (77) 



Note that Eq. ( 75 1 can be inverted to give 



a= i=2 {bo 



2 i(^- 6 °)+\/V 61 ' 



^=y^^- b ^-y^- ^ 

From the Hamiltonian, we easily find the time evolution of the operators 
bo, b\ and b 2 in the Heisenberg picture, i.e. 

60 (t) = e-^- Q n (Q), 
b 1 (t)=e- i » t b 1 (0), 

b 2 (t) = e- l{uJ+n)t b 2 (0). (79) 
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Combining these relations with ( 75 1 and ( 78 ) , we obtain the time evolution of 
the field and atomic operators: 



a(t)-- 



cos(fM)a(0) — i 
k 



sin(Oi)^i(0) 



/ _/V _ ^ 
/\/^^sin(m)6(0) 



y sin(fit)a(0) + (-^ r 



7V 



AT 

3(m))a(o) 



TV 



s(fit) - 1)6(0) 



N 
N-k 
N 



k , n , , , Jk(N-k) 



iV 



cos 



(«*))6(0) 



3 (nt) - 1)^(0) 

(80) 



Let us first consider the entanglement condition in Eq. ([2| , and choose A = £j 
and B = ^2- The resulting inequality will tell us if the two groups of atoms are 
entangled. If the initial state is of the form 



|*> = M|0,0> ! 



(81) 



then the quantities at time t appearing in the entanglement condition can be 
easily calculated, since only the terms involving the field operator at t — can 
be non-zero. In particular, we obtain 



(&(*)&(*)) 
(£(*)&(*)) 



N 



sm z (Qt) (i>\(a(0)y 



k 

N 
N 



sin 2 (m) (^1^(0)0(0)1^), 
~ fc sin 2 (m) (^| O t(0)o(0)| 



N 



hence for t > 0, the entanglement condition becomes 

|(V»|(a(0)) a |^)| 3 >(^|ot(0) O (Q)^) a . 



(82) 



(83) 



We first note that a coherent state will not satisfy this condition. That is not 
surprising, because under the action of our Hamiltonian, an initial state with 
the field in a coherent state and the atoms in their ground states will evolve 
into a product of coherent states in all three modes. Such a state is clearly not 
entangled. If the field is initially in a squeezed vacuum state = S(z) |0) the 
situation is very different. The entanglement condition reduces to 



cosh 2 r > sinh 2 r 



(84) 



which is satisfied for any r ^ 0. Therefore, squeezing in the initial field will lead 
to entanglement between the two groups of atoms. 

Let us now consider the entanglement condition in Eq ([TJ and again choose 
A = £i and B = 6- If the initial state at t — is of the same form as before, 
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\^f) = 1 0, 0) , then the condition can be again greatly simplified as most of 
the terms on both sides of the equation are zero. We find that 




(fit) (^|(at(0)a(0)) 2 |^) 



(fit) (V|at(0)a(0)M, 



(85) 



and the entanglement condition becomes 



(^| a t( )a(0) |^) 2 > (at(O)a(O)) 2 |V> - <VI a f (0)o(0) |V) ■ 



(86) 



Expressed in terms of the photon number operator, n = a^a, this is just 
A 2 n(0) < (n(0)), where A 2 n = (n 2 ) — (n) 2 is the variance of n. Therefore, 
we see that there will be entanglement between the two groups of atoms pro- 
vided that the field is initially in a state with sub-Poissonian photon statistics. 
Connections between nonclassical states and entanglement have been noted be- 
fore. For example, for the output state of a beam splitter to be entangled, the 
input state must be nonclassical [12 . The amount of two-mode entanglement 
that can be produced by a single-mode field using linear optics, auxiliary clas- 
sical fields and ideal photodetectors has even been proposed as a measure of 
the nonclassicality of a state [13]. Here we see yet another connection between 
nonclassical states an entanglement in the dynamics of the Dicke model. 

Light has, in fact, been used to entangle atomic ensembles, though this 
is often done by letting the light interact with both sets of atoms and then 
measuring it []3]-[T7]. The system we have considered here is a cartoon version 
of the one in the experiment |17j . but it does provide useful information in that 
it shows that certain kinds of nonclassical light can lead to atomic entanglement. 

6.2 Two beam splitters 

Let us now consider a system of field modes that is closely related, at least in 
its description, to the atom-field system we have just studied. The system is 
depicted in Fig. 2. It consists of three field modes and two beam splitters. One 
mode, mode a, passes through both beam splitters, while the other two, modes 
b and c, pass through just one. Our goal is to entangle modes b and c by making 
them interact with mode a through the array of two beam splitters. Because 
this system closely resembles our previous example where a light mode was 
used to entangle two groups of atoms, it is not surprising that the entanglement 
condition Eq. (JTJ again provides a usefull tool to derive simple conditions on 
the input state of the ancillary mode a that will result in modes b and c being 
entangled at the output. This condition was applied to determine when the 
output of a single beam splitter is entangled in |18j . 




The system we are examining here was first studied in |19j . In that paper a 
series of beam splitters was used to model a reservoir, and, conditions for the 
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entanglement of reservoir modes, which would correspond to our modes b and 
c, were found. However, in that paper only Gaussian states were considered. 
We recall that a beam splitter acts on two input modes a and b as 



dout 

bout 



tai n + rh 
—rat 



in ■ 

tbjn 



(87) 



where r and t are positive, and r 2 +t 2 = 1. For beam splitters with transmittance 
and reflectance (ii,J"i) and (^2,^2) respectively, the relationship between the 
output modes and the input modes is given by 

dout 



C-out 



-riitidin + nhn) + t 2 Ci, 



(88) 



Let us suppose for simplicity that modes b and c are initially in the vacuum, 
i.e. l^in) = \ip) |0) |0). The entanglement condition in Eq. ([I]) for A = b out and 
B = c ou t is 

> ■ (89) 



Replacing b out and c ou t by their expression Eq. (|88| we find 

tirir 2 



aLa,i 



Win a in) \$) 



! aL at- 



{°out c out) 

(blut b outcl ut c out ) = (r 1 t 1 r 2 )' 2 
and the entanglement condition becomes, assuming that r\t\ri ^ 0, 
{{d\ n d in f) - {a\ n a m ) 2 < {a\ n d in ). 



W) > (90) 



(91) 



i.e., b and c will be entangled provided that mode a is initially in a state with 
sub-Poissonian statistics. 
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Now let us see what happens when we apply a more powerful entanglement 
condition. Motivated by considerations similar to those in Section V, we set 



1|. We then find that a state is 



A = ziKut + z 2b ut and B = c out in Eq. 
entangled if the 2x2 matrix M, where 

I — (boutb ut c out c out) > 

" out C out C out) 



Ma 
M 22 



\{b out i 

M 21 = (boutC ut)(boutct ut ) - (b 2 Q •■' 



= \{b\ ut Cout)\ 2 - (btutbout 



,t 



utC'out c out) , (92) 

has a positive eigenvalue, which will be the case if |Mi2| 2 > M\\M 22 . When we 
express this condition in terms of the input operators, and assume that the b 
and c modes are initially in the vacuum state, we find that the two- mode output 
state is entangled if 



\(al)(n)-(nal)\ 2 >[(n)-(A 2 n)] 



\(al) 



1 - 



1 



(93) 



where n — 



Note that the Schwarz inequality implies that the second 



factor on the right-hand side is negative, so that if the input field has sub- 
Poissonian statistics, then this inequality is satisfied. Therefore, this condition 
includes the one in Eq. (91 1. There are, however, states that satisfy the new 



condition but do not satisfy the old one. For example, for the state 

1/2 



1 



1 

vi 



i 

71 



|2>, 



(94) 



we find that (A 2 n) — (n) = — 2v2e, to lowest order in e. On the other hand, 



Eq. (93 1 becomes, again keeping only lowest order terms in e 



1 



> 



-V2e. 



(95) 



So, if e < and small, then the state does not have sub-Poissonian statistics, 



but Eq. (93 1 is satisfied, so the output will be entangled. Therefore, Eq. (93) is 



a stronger condition than Eq. (91 ) 



7 Local uncertainty relations and two modes 

So far we have been concentrating on the consequences of the inequalities in 
Eqs. ^ and Needless to say, it is possible to derive other entanglement 
conditions involving non-hermitian operators. One way, which we will briefly 
explore here is to apply a variant of the local uncertainty entanglement condition 
due to Hofmann and Takeuchi [20 . Let us start by considering two modes. Now 
suppose that A is an operator on mode a and B is an operator on mode b. We 
want to first find an expression for ((A^ +B jf )(A+B)) — \ (A+B) \ 2 for a separable 
state 

P = T.P^ ] ®Pi b) - (96) 
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We first note that for an operator D and a density matrix 

P=^2pkPk, (97) 

k 

we have that 

k 

> Y, Pk ((rtD) k -(rt) k {D) k ), (98) 

k 

where expectation values with the subscript k denote the expectation value with 
respect to p k . If we now apply this to D = A + B with the separable density 
matrix above, we find that 

((A^ + B^(A + B))-\(A + B)\ 2 >Y,Pk((^A) k -\(A) k \ 2 + (B^B) k -\(B) k \ 2 ). 

k 

(99) 

This condition can easily be extended to the case in which we have more than 
one operator for each subsystem. If we have M operators for mode a, Aj and 
M for mode b, Bj, then the above inequality generalizes to 

M 



+ B]){A, + Bj)) \{Aj + Bj)f > 

M 

Y^Pk Y,({A]Aj) k - \{A,) k \ 2 + {B)Bj) k \{B,) k \ 2 ). (100) 
k j=i 

These two inequalities will hold for all separable density matrices, so if a state 
violates them, it must be entangled. 

As a simple example of an entanglement condition that can be derived from 
Eq. (99 1, let us set A = a and B = tf . We then have that 

(( t + & )( a + 6 t))_|( + 6 t)|3 > ^j, fc ((oto) fc -|(o) fc | 2 + (66t) fc -|(6) fc | a ) 

k 

> 1- (101) 

If this inequality is violated by a state, that state is entangled. A two-mode 
squeezed vacuum state, 

I^He^-^lO), (102) 
will violate this inequality. We find that for this state 

((^ + b)(a + &t)) - |( a + 6 f )| 2 = e~ 2r , (103) 

so that for r > 0, the condition is violated, and we can conclude that the state 
is entangled. It should be noted that this condition can also be derived from 
the analysis given by Shchukin and Vogel [8]. 
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A very similar condition can be derived for atom- field entanglement, 
have N two-level atoms represented by collective spin operators S + , S~ 
S z , then setting A = and B = J + we find that for a separable state 



If we 

, and 



((a + S-)(a^ + S+)) -\(a + S-)\ 2 > 1, 



(104) 



so that any state that violates this inequality must be entangled. A very simple 
example is the state, for N = 1 



for < 9 < tt/4. 

8 Conclusion 

We have presented a number of entanglement conditions that employ non- 
hcrmitian operators. While these conditions follow from the partial transpose 
condition, they are much simpler to use. It is only necessary to compute a small 
number of correlation functions rather than to have to diagonalize the partially 
transposed density matrix for the entire system. We used these conditions to 
show the presence of entanglement in a number of systems of interest in quantum 
optics, including Jaynes-Cummings model and the Dicke model. The conditions 
are sufficiently flexible that we can study entanglement between modes, entan- 
glement between a field mode and an atom, and entanglement between groups 
of atoms. Most of these conditions are invariant under some set of local unitary 
transformations. For field modes we considered local Gaussian operations. By 
making the conditions invariant under a set of local operations, we are able to 
make them able to detect entanglement in a wider class of states. 
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Appendix 

Here we present the details needed for the calculations involving the Tavis- 
Cummings model. As we noted, the Hilbert space splits up into a direct sum 
of invariant subspaccs, , which are characterized by the total excitation 

number, n. is spanned by {\0,g,g)} so the ground state of the system is 



|V>) = cos0|e)|O) +sin%)|l) 



(105) 




(106) 
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with eigenvalue 



4 0) = (0,g,g\n\0,g,g) = -LO. 



(107) 

is spanned by {\l,g,g) , \0,e,g) , \0,g, e)}, and the Hamiltonian can be 
written in this basis as 

k k \ 

= | k . (108) 

K / 

This Hamiltonian has three eigenvalues: 



(1) 



(109) 



to which we associate the three eigenvectors 



1>\ 



(1) 



V2 



[|0, e, ff > - |0, ff ,e>] s 



(110) 



is spanned by {|n, g, g) , \n — 1, e, 3) , |n — 1, e) , |n — 2, e, e)}, and the 
Hamiltonian can be written in this basis as 



/ w(n-l) 







V 



K-^n oj(n — 1) fiVn — 1 
K-^/n u>(n — 1) — 1 

K\/n — 1 K\/n — 1 cj(n — 1) J 



(111) 



This Hamiltonian has four eigenvalues: 



E. 



,(n) 


,(n) 
1,2 

(n) 



1) - Kv /2(2^T), 
Lo(n — 1), 

1) + Ky/2(2n-l), 



(112) 



to which we can associate the four eigenvectors 



(n) 



n 1 r |n- l,e,#) + |n - l,.g,e)- 
2(2^1) I"' ff> - ^ [ 71 



/ n- 1 
2(2ra- 1) 



|n — 2, e, e) 



(n) 



>/2 



[|n- l,e,.g) - \n - 1, g, e) ] . 
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4 n) 



4 n) 



n- 1 
2n- 1 



|n-2, e,e), 



2(2n- 1) 

/ n-1 
V 2(2n - 1) 



2ra- 1 

1 r|n-l,e,g) + |n - 1,5, e) 
v^ 1 



V2 



|n-2,e,e) 



(113) 



Suppose that at time t = 0, the field contains n excitations while both atoms 
are in their ground state, i.e. 



|*(0)) = \n,g,g) 



(114) 



2(2n- 1) 



(n) 



n-1 
2n - 1 



4 n) 



2(2n- 1) 



Note that this expression is valid for any n > 0, provided that the appropriate 



eigenvectors are chosen. In particular, when n = 1, 



as defined in Eq. ( 110 1 



«/>3 ) corresponds to 



4» 



After a time i, the state of the system will have evolved according to 



!*(*)) 



*(0)> 



iE K ' >t 



2(2n- 1) 



e '"o 



(») 



n — 1 
2n - 1 ' 



iE\ >t 



-e 



(n) 



2(2n - 1) 



-2u(n — 1)£ 



l 



— (ncos(fii) + n-l) |n,5,fli> 
-2n — 1 

^n(n - 1) 



(2n - 1) 



cos(fH) - 1) |n-2,e,e) 



(115) 



where il = Ky / 2(2n — 1). 



References 

[1] R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Rev. Mod. 
Phys. 81, 865 (2009). 

[2] O. Glihne and G. Toth, Physics Reports 474, 1 (2009). 

[3] R. Simon, Phys. Rev. Lett. 84, 2726 (2000). 

[4] L. -M. Duan, G. Giedke, J. I. Cirac, and P. Zollcr, Phys. Rev. Lett. 84, 
2722 (2000). 



23 



[5] S. Mancini, V. Giovannetti, D. Vitali, and P. Tombesi, Phys. Rev. Lett. 
88, 120401 (2002). 

[6] G. S. Agarwal and A. Biswas, New Journal of Physics 7, 211 (2005). 

[7] Mark Hillcry and M. Suhail Zubairy, Phys. Rev. Lett. 96, 050503 (2006). 

[8] E. Shchukin and W. Vogel, Phys. Rev. Lett. 95, 230502 (2005). 

[9] Geza Toth, Christoph Simon, and Juan Ignacio Cirac, Phys. Rev. A 68, 
062310 (2003). 

[10] Hyunchul Nha and Jacwan Kim, Phys. Rev. A 74, 012317 (2006). 

[11] V. Buzck, M. Orszag, and M. Rosko, Phys. Rev. Lett. 94, 163601 (2005). 

[12] M. S. Kim, W. Son, V. Buzck, and P. L. Knight, Phys. Rev. A 65 032323 
(2002). 

[13] J. K. Asboth, J. Casamiglia, and H. Ritsch, Phys. Rev. Lett. 94, 173602 
(2005). 

[14] L. -M. Duan, J. I. Cirac, P. Zoller, and E. S. Polzik, Phys. Rev. Lett. 85, 
5643 (2000). 

[15] A. Kuzmich and E. Polzik, Phys. Rev. Lett. 85, 5639 (2000). 

[16] For a review see K. Hammerer, A. S. Sorensen, and E. S. Polzik, quant- 
ph/0807.3358. 

[17] B. Julsgaard, A. Kozhckin, and E. Polzik, Nature 413, 6854 (2001). 

[18] Mark Hillcry and M. Suhail Zubairy, Phys. Rev. A 74, 032333 (2006). 

[19] Daniel Nagaj, Peter Stelmachovic, Vladimir Buzek, and Myungshik Kim, 
Phys. Rev. A 66, 062307 (2002). 

[20] Holger Hofmann and Shigcki Takcuchi, Phys. Rev. A 68, 032103 (2003). 



24 



